Gonihedric (and Fuki-Nuke) Order 



D. A. Johnston 

Dept. of Mathematics and the Maxwell Institute for Mathematical Sciences, 
Heriot-Watt University, Riccarton, Edinburgh, EH14 4AS, Scotland 

Abstract. A 3D Ising model with a purely plaquette, 4-spin interaction displays a 
planar flip symmetry intermediate between a global and a gauge symmetry and as a 
consequence has a highly degenerate low temperature phase and no standard magnetic 
order parameter. This plaquette Hamiltonian is a particular case of a family of 3-D 
Gonihedric Ising models defined by Savvidy and Wegner. An anisotropic variant of the 
the purely plaquette Gonihedric model, originally discussed as the "Fuki-Nuke" model 
by Suzuki, is non-trivially equivalent to a stack of 2D Ising models, each of which can 
magnetize independently at the phase transition point. 

Consideration of this anisotropic model suggests that a suitable order parameter 
in the isotropic case may also be constructed using a form of planar magnetization, 
in which nearest neighbour correlators {&i<Jj) summed over planes replace a sum over 
spin values (<7i). We conduct Monte-Carlo simulations to investigate this and related 
candidate order parameters in a toy model of Gonihedric ground states, the Fuki-Nuke 
model and the isotropic plaquette Gonihedric model itself. 
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1. Introduction 

The Gonihedric Ising model was defined by Savvidy and Wegner [I] as a cubic lattice 
transcription of an earlier Gonihedric random surface model developed by Savvidy [2] 
which was intended to be a discretization of string theory. In the case of the Gonihedric 
Ising model, the spin cluster boundaries were used to model a gas of surfaces and the 
spin couplings were chosen so that edges and intersections of plaquettes rather than 
their surface areas were weighted, which was also the defining feature of the Gonihedric 
random surface model. 

Such a weighting can be arranged using a one-parameter family of generalized 3D 
Ising models with fine-tuned nearest neighbour (ij), next to nearest neighbour ((ij)) 
and plaquette interactions [ijkl]. The Hamiltonians in this family are given by 

Hgonihedric = — 4/C ^ CTiOj + K ^ O \0 j — (1 — ft) OiOjO k Ol , (1) 

(ij) ((ij)) [ijkl] 

where the parameter k gives the relative weight of two plaquettes meeting at right 
angles and four-plaquette intersections in the spin cluster boundaries. One member of 
this family, the purely plaquette action with k = 

H K= o = - OjOj(y k oi (2) 

[ijki] 

has merited particular attention due to its interesting properties, both static and 
dynamical (3} [3]. It has a strong first order transition surrounded by a region of 
metastability and shows evidence of glass-like behaviour at lower temperatures in spite 
of having no quenched disorder. 

The symmetry properties of both the ground states and low temperature phase of 
H K=0 are unusual. The standard Ising model with only nearest neighbour interactions 

H ising = — Q- j0~ j (3) 

(ij) 

displays a twofold symmetry in its ground state, whereas a 3D Z 2 gauge theory (where 
the spins live on plaquette edges rather than vertices as in H K=0 ) 

Hcauge = ~ UijUjkUklUu (4) 

[ijki] 

has the local gauge symmetry — >■ SiUijSj, with Sj = ±1. The symmetries of plaquette 
Gonihedric model of equ. ^ lie somewhere between these since it is possible to flip 
arbitrary, possibly intersecting, planes of spins in the ground state at zero energy cost. 
Both Monte-Carlo simulations and low temperature expansions [5] confirm that this 
remains a symmetry of the low temperature phase, which is thus highly degenerate. 
A consequence of this symmetry is that the standard magnetization is zero in the low 
temperature phase. The nature of the magnetic order has thus remained unclear. 

In this paper we discuss the recent suggestion by Hashizume and Suzuki [6] that a 
suitable order parameter for H K=0 might be constructed using nearest-neighbour spin- 
spin correlations (o~iO~j), which was prompted by Suzuki's earlier study of an anisotropic 
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Gonihedric model [7]. We conduct measurements of this candidate order parameter on 
a toy ensemble of ground state configurations constructed by flipping planes of spins 
starting from a ferromagnetic starting configuration to mimic the ground state structure 
of the Gonihedric model, before moving on to simulations of both the anisotropic and 
isotropic Gonihedric models themselves. 



2. The Anisotropic Gonihedric Model (a.k.a. the Fuki-Nuke Model) 

The study of a strongly anisotropic variant of the plaquette Gonihedric Hamiltonian by 
Suzuki [7] pre-dates the work of Savvidy and Wegner by some twenty years. For reasons 
that will become apparent, Suzuki dubbed this the "Fuki-Nuke" or "no roof" model 
and showed that it was (non-trivially) equivalent to a stack of 2D Ising models. This 
was later rediscovered by Jonsson and Savvidy [8] using a transfer matrix approach and 
discussed by Castelnovo et.al. (9] using methods more akin to Suzuki's original approach. 

In Suzuki's work the principal concern was to obtain solvable 3D spin models, and 
the Fuki-Nuke model was constructed as a specific example. The starting point was an 
anisotropically coupled variant of H k= q 

H an iso = -Jx 2j ViVjVkVl ~ Jy / ] ^i^jCTk^l ~ J Z / ] Cri^jOk^l (5) 
[ijkl] [ijkl] [ijkl] 

where J x ,J y , Jz were the couplings for plaquettes perpendicular to the x , y , z axes 
respectively. Setting J x = J y = 1 and J z = gave a model where the horizontal 
"roof" plaquettes had zero coupling, hence the sobriquet of "Fuki-Nuke". 

To show that this is equivalent to a stack of 2D Ising models requires a little further 
work, which can be done succinctly by using the non-linear er-r transformation, following 
Suzuki [7] and Castelnovo et.al. [9j. We first define bond spin variables Tj = OiOj that 
can be thought of as living at the centre of (only) the vertical lattice bonds. Since we 
have condensed the site labels on both sets of spins, we write a k and r k to represent any 
of the spins in the kth horizontal 2D layer. With this notation the direct transformation 
for a spin in the kth layer is 

Ok = TiT 2 n ■ ■ ■ T k (6) 

and the inverse transformation is given by 

Ti = 0i, t 2 = (Jia 2 , . . . , t n = a N -i<j N (7) 

In order to obtain a one-to-one correspondence between a and r spin configurations we 
need to specify the value of the <j\ / T\ spins on a given (first in this case) horizontal 
plane, but in all the other planes a product of two a spins at the end of each vertical 
bond gives the r spin. The resulting Hamiltonian is then 

H = - { T 'iTi+e x + nTt+eJ (8) 

i 

where e x , e y are unit vectors in the horizontal directions. We have not written the 
contribution of the fixed plane, which will vanish in the thermodynamic limit. The T{ 
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spins are coupled to their horizontal nearest neighbours in the e x and e y directions and 
there are no vertical inter-plane couplings. We have thus shown that the Fuki-Nuke 
model is equivalent to a stack of uncoupled 2D Ising models arranged horizontally one 
above each other. 

The magnetic order parameter for the ith 2D Ising layer is given by the standard 
expression 



M ^ = (jf2 £ T ) ( 9 ) 

\ single plane / 

which may be rewritten in terms of the original Oi spins as 



M 2D)i =l— 2 £ ff '^) • ( 10 ) 

\ single plane I 

where e z is the unit vector in the vertical direction. The nearest neighbour <7j<7j + e z 
correlator is measured on all the vertical bonds bisected by the horizontal plane of the 
corresponding n spins. The unusual appearance of a two-spin correlator as an order 
parameter is simply the result of expressing the magnetization in terms of the original 
oi variables rather than the bond spins 7$. Remembering that T\ = <Ji<Ji + e z and using 
the standard results for the 2D Ising model magnetization we have [7] 

M 2D ,i = £ = ± (1 ~ smh- 4 (/3)) 1/8 (11) 

\ single plane I 

which will behave as ± \/3 — (3 C \ 5 near the critical point f3 c = ln(l + y/2). More generally, 
since = T\ r 2 . . . and the different n layers are decoupled, 

M 2D ,i,n = (jjj £ (?i(Ti + ne z J = {M 2D)i ) n . (12) 
\ single plane I 

The in-plane horizontal correlations between "dimers" of <7j spins, where the i,j spins 
now lie on the same plane, 

((Vi Vi+eJicTj a j+ e z )) = fa Tj) (13) 

are also, of course, Ising- like. 

Since the various 2D Ising layers magnetize independently we could construct a 
quasi-3-D order parameter to be of the form 

k £ (\w S a ^+A) (14) 



M abs 

xy planes single plane 



or 



xy planes \ single plane 
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to avoid inter-plane cancellations. We have explicitly retained the various normalizing 



factors in equs.(14, 15) for a cubic lattice with iV 3 sites. We could also define similar 



quantities for non-nearest neighbour correlators 

M abStn = ^ E (|^2 E Wi+m,]) (16) 

xy planes single plane 

or 

M sq>n = ± E ((^ E ) 

xy planes \ single plane / 

with n > 1. Since each plane magnetizes at the same point the resulting expressions are 
given by 

sinh- 4 (/3)) 1/8 

„-4/m\l/4 



M abs 


= (i- 


M sq 


= (i- 


M^abs, n 


= (i- 


M-sq, n 


= (i- 


i other 


variants 



sinh- 4 (/3))" /8 
sinh- 4 (/3)) n/4 . 

the magnetization are given in terms of M abs = \M2D,i\ 
this is probably the most suitable choice of "three-dimensional" order parameter for 
the Fuki-Nuke model. In addition, it emphasizes that the critical behaviour of the 
Fuki-Nuke model is really that of the 2D Ising model in each layer. 

In the anisotropic model the orientation of the planes in the summation is fixed 
by the choice of zero coupling, which we have taken to be horizontal planes in our 
case. The suggestion by Hashizume and Suzuki in [6] is that a similar order parameter 
measuring coplanar, nearest neighbour spin correlations could still serve for the isotropic 
Gonihedric model. They considered a mean-field approach to the model and correlation 
function inequalities in order to corroborate this and found supporting evidence from 
both, as well as extending the discussion to disordered couplings. Here, we investigate 
the order parameter in the isotropic 3D Gonihedric model directly using Monte-Carlo 
simulations, restricting discussion to the purely ferromagnetic case for simplicity. 

In equs. ( 14|15 ) we have introduced either a modulus or square for the spin 



correlators in order to avoid inter-plane cancellations, whereas Suzuki's approach is 
tantamount to considering the direct equivalent of equ. ^ 

M = ^E(^ E *<*h 

planes single plane 

in the presence of an external field (later taken to zero) which picks out one of the 
2D magnetizations on each plane. The possibility of using an order parameter akin to 



the M sq in equ. ( 15 ) has been mentioned previously by Lipowski [3] , who noted in a 
simulation of H k= q that it appeared to possess the correct behaviour. He also remarked 
on the possibility of using non-nearest neighbour correlators as the order parameter (s), 
which we shall discuss further below. In the remainder of the paper we compare 
measurements of the candidate order parameters on toy ground state configurations, 
the Fuki-Nuke model and the isotropic 3D Gonihedric model itself. 
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3. Monte Carlo investigations: Toy Model, Fuki-Nuke and the Real Thing 

3.1. Toy Model - Ground States 

Ground states of the Gonihedric Hamiltonian H K=0 may be obtained by flipping some 
set of planes of spins in a reference ferromagnetic configuration which has all spins up or 
down. A simple test for any candidate order parameter is then to generate an ensemble of 
such configurations and measure its expectation value. A typical configuration on a finite 
cubic lattice generated in this manner would look something like that in Fig. (1). We can 




Figure 1. A "ground-state" configuration generated by flipping arbitrary planes of 
spins, shown dotted, starting in a ferromagnetic state with all the spins positive or 
negative. 



also measure proposed order parameters on purely random (i.e. infinite temperature) 
configurations. Any that fail to distinguish between the ground states and random 
configurations would not, at first sight, appear to be suitable for use in a simulation of 
the full Hamiltonian. 

If the suggestion of ||6j is correct, we would expect that the order in the low 
temperature phase of the isotropic 3D Gonihedric model could be discerned using an 
appropriate generalization of either the M a b s and M sq in equs. ( I4p5 ). For isotropic 
configurations such as that in Fig. (1) the system should be agnostic to the orientation 
of the planes used in calculating the correlators, so any of the three choices of orientation 
for reference planes in 

Yl a i a i+e x , y , z ) ( 19 ) 



M xy,z 

abs 



1 

N 



E 

yz\xz\xy planes 



N 2 



single plane 



or 



M x,y,z 



1 

N 



E 

yz\xz\xy planes 




single plane 



(20) 



where e x , e y , e z are the unit vectors in the x, y, z directions respectively, should be able to 
serve as order parameters. In all cases the correlators are nearest neighbour correlators 
for the (jj spins on bonds perpendicularly bisected by the reference planes, as in the 
Fuki-Nuke model. Reassuringly, we find that the various M^ y s ,z from equ. ( 19 ) and the 
M*' y,z from equ. (20) all give +1 when measured on an ensemble of flipped states such 
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as that in Fig. (1) so they do, indeed, pick up the planar order that is characteristic of 
the Gonihedric models. In addition, the M^' z and M^ y ' z all give zero when measured 
on purely random configurations, so they are capable of distinguishing this order from 
the complete disorder at infinite temperature. 

The choice of nearest neighbour correlators in defining the order parameter does 
not appear to be obligatory, though given the line of reasoning leading to it from the 
Fuki-Nuke model it is probably the most natural. If we measure quantities such as 

^1,11 = ^ E E Wi + ne x , y , z \) (21) 

yz\xz\xy planes single plane 

or 



2 



M stn E ( \W2 E ) ( 22 ) 



N 

yz\xz\xy planes \ single plane 

with n > 1 on the ensemble of flipped ground state configurations these also give +1 
(and zero on random configurations). That this should be so is clear from starting with 
configurations in which planes of one orientation only, say horizontal, are flipped. For 
these, plane-to-plane correlations at any separation must be ±1. Taking a modulus or 
square for each planar sum then makes all of these +1. Flipping planes of spins in the 
two remaining directions perpendicular to this will change the signs of lines of spins on 
both the horizontal planes contributing to the correlator sums and therefore not change 
the values of the contributing correlators. 

When the standard magnetization is zero 

M = ^E^) = ° (23) 

i 

the magnetic susceptibility x is given by 

iJs E E^> = ^ E E^+«> (24) 

i j in 

which is related to the class of order parameters discussed here, since the sums may be 
arranged in a similar manner to those in M sq and M abs . However, the expression for 
X includes a summation over distances greater than one and also retains the signs of 
the various different contributions in the sum. Consequently, the susceptibility does not 
distinguish between configurations such as those in Fig. (1) and random configurations, 
giving +1 in both cases, but as we shall see below it does appear to display the 
characteristics of an order parameter in a simulation of the full isotropic Gonihedric 
model. 



X N 3 



3.2. Fuki-Nuke Model 

Since the anisotropic Gonihedric Hamiltonian, or Fuki-Nuke model, originally inspired 
the definition of the candidate order parameters it too can serve as a test case for 
simulations. We know that 2D Ising criticality is expected in each horizontal layer, 
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so the transition point will be seen in the continuum limit at the 2D Ising value of 
(3 C = ln(l + a/2) ~ 0.88. In the Fuki-Nuke model the orientation of the the 2D planes 
used to measure the nearest neighbour spin correlations is fixed by the choice of the 
zero couplings. In our simulations we take J z = 0, which gives horizontal planes in the 
sums and means that M z bs should be the correct choice for the order parameter. 

We simulate the Fuki-Nuke model on a 20 3 lattice with periodic boundary conditions 
using 10 7 Metropolis Monte-Carlo measurement sweeps with a hot start and an 
equilibration time of 10 5 sweeps as a test case. In Fig. (2) we can see that both M a b s 
and M sq give a clear signal for the pseudo-critical point, which can be observed close to 
the continuum value of ~ 0.88 already on the 20 3 lattice via, for example, the peak in 
the specific heat measurements. As a consequence of cancellations between differently 
signed 2D layers of Ising spins in the magnetized phase the standard magnetization M 
remains zero for all in the Fuki-Nuke model and does not give an obvious signal for 
the transition. 




0.2 0.4 0.6 0.8 1 1.2 



Figure 2. The energy and order parameters M a }, s and M sq on a 20 3 lattice for the 
Fuki-Nuke model. 



3.3. (Isotropic) Gonihedric Model 



The preceding measurements reported for the toy ground state model and the Fuki-Nuke 



model suggest that any of M x ab f from equ. (19), the M^ y,z from equ. (20) or even the 



^■abs'n an d ^sq,n °f equs. ( 21~|22 ) might be viable candidates for an order parameter in 
the isotropic 3D Gonihedric model itself. To clarify this, we carried out simulations of 
the isotropic Gonihedric model on 10 3 and 15 3 lattices using 10 7 Metropolis Monte-Carlo 
measurement sweeps, following a hot start and an equilibration time of 10 5 sweeps. The 
various M^ y s ' z and M^ y ' z along with M x ab f n and Mf* z for n > 1 were measured on 
both lattices. 

In Fig. (3) we plot the measurements of M^ y,z and M*l v ' z on the 10 3 lattice 



abs 



and in Fig. (4) on the 15 lattice. The measurements of the M's for the different 
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x, y, z orientations of the reference planes are indistinguishable from each other and 
the observed jump in the various M's coincides with the pseudo-critical values of 
(3 determined from the jump in the energy observed in simulations (~ 0.54 on the 
10 3 lattice and ~ 0.57 on the 15 3 lattice with "hot" starts). These results confirm 
that the orientation of the reference planes is, indeed, irrelevant in measuring these 
magnetizations and that all three choices are equivalent. They also show clearly that 
the correlators distinguish between the disordered and ordered phases of the isotropic 
model as suggested by the arguments in 



0.9 - 
0.8 - 
0.7 - 
0.6 - 
2 0.5 - 
0.4 - 
0.3 - 
0.2 
0.1 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

P 

Figure 3. M^' z (upper) and M^ v ' z (lower) on a 10 3 lattice for the isotropic 
Gonihedric model. Lines are drawn to guide the eye. 
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0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 

P 

Figure 4. Mabs* ( u PP er ) an d M^ y,z (lower) on a 15 3 lattice for the isotropic 
Gonihedric model. Again, lines are drawn to guide the eye. 

To guard against the possibility of discarding too much information by considering 
a simple global order parameter we can measure the distribution P(M 2 D,i) of the 
coplanar 2D spin correlations. In Figs. (5,6) we have combined the three orientations 
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of reference planes to plot a P(M2d) in the high temperature and low temperature 
phases respectively. The three individual orientations contributing to this give identical 
distributions. In Fig. (5) P(M2d) is gaussian (narrowing with system size and increasing 
temperature) and centred at M 2 d = 0, whereas in the low temperature phase in 
Fig. (6) at j3 — 0.58 there are two sharp peaks at M 2 d = ±1- We would expect the 
slight asymmetry between the peaks in the plot to disappear with better statistics and 
increasing lattice size. The observed behaviour strongly supports the hypothesis of Fuki- 
Nuke type order, where coplanar correlators "magnetize" independently of the other 
planes, in the isotropic model too. The plots also demonstrate that a [0, 1] type order 
parameter is sufficient for the isotropic Gonihedric model, since there is no non-trivial 
structure in P(M 2 d) such as that seen in the overlap distribution P(q) of (mean-field) 
spin glasses. 

0.7 
0.6 
0.5 
0.4 

0.3 
02 
0.1 


-0.6 -0.4 -0.2 0.2 04 0.6 

u 

Figure 5. P(M 2 d) at /3 = 0.20 in the high temperature phase on a 10 3 lattice 




3.5 
3 

2.5 



1.5 




Figure 6. P(M 2 d) at /3 = 0.58 in the low temperature phase on a 10 3 lattice 

We also find, as with the toy ground state ensemble, that M^ y s ' z n and M^rf for 
n > 1 show very similar behaviour to M^f z and M^' y,z . In Fig. (7) M* bs 6 on a 15 3 
lattice is plotted, along with M z ahs for comparison. M* bs 6 is typical of the various M a b StTl 
(and M sq>n ) for other n. As might be expected the magnetization constructed from 
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Figure 7. M* bs (upper) and M* bs 6 (lower) on a 15 3 lattice for the isotropic Gonihedric 
model. 



the more widely spaced correlators in M^ fes 6 is lower in the disordered phase, but it too 
jumps sharply to one at the transition point just as does M* bs . Since we employ periodic 
boundary conditions the measured magnetizations repeat from M a b s j on the 15 3 lattice 
as the distances in the correlators wrap around the lattice. With the isotropic model 
there is no longer a mapping to the 2D Ising model, unlike the Fuki-Nuke model, so we 
do not have an exact expression for the planar average of nearest neighbour correlators 
as in equ. ( 11 ), nor do we have any assurance that different (bi-)layers of spins decouple, 
which allows us to express everything in terms of the nearest neighbour correlators 
alone in the Fuki-Nuke case. Nonetheless, it would appear that the most fundamental 
correlations in the isotropic case are also those between nearest neighbour spins. The 
observed signal for the transition in the higher correlators are a consequence of that in 
the nearest neighbour correlations, just as with the toy ground states and the Fuki-Nuke 
model. 

A curious feature of the isotropic Gonihedric model simulations that should be 
highlighted is that the standard magnetic susceptibility, as defined in equ. (24): x — 
-^3j2 i J2j(< T i< T j) = F Ei En^+n)' a ^ so behaves as an effective order parameter, 
unlike the toy ground state ensemble. This can be seen in Fig. (8) where \ is +1 at 
high temperatures and drops sharply to zero at the observed pseudo-critical points. 
We have already seen in our discussion of the toy model for ground states that \ does 
not distinguish between random, infinite temperature configurations and the flipped 
ground-state configurations, which means that the observed behaviour is not a direct 
consequence of \ detecting the coplanar order of the low temperature phase as with 
^abslq- ft appears rather to be a dynamical ergodicity breaking effect. At high 
temperatures the flip symmetry of the Hamiltonian is not broken, so \ must be one, 
whereas the observed value of zero for \ throughout the low temperature (high 0) 
phase is due to the system "freezing" into one configuration as the Metropolis spin flip 
acceptance drops rapidly to zero at (3 C . It would be interesting to explore the behaviour 
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Figure 8. The magnetic susceptibility \ 011 both the 10 3 and 15 3 lattices is one at 
high temperatures and drops sharply to zero at the observed pseudo-critical points. 



of x in a multicanonical simulation such as that carried out in 10 . In this case the 
system would still be free to explore multiple ground states, so the behaviour of x would 
presumably be more akin to that in the toy ground state ensemble (and less like an order 
parameter). 

It is also worth noting that the strong hysteresis associated with the first order 
transition means that a more accurate estimate of the continuum transition point /3 c (oo) 
is best done with such a multi-canonical simulation, or similar. In [10] the value of 
/3 c (oo) = 0.54925(6) obtained from the finite size scaling of the specific heat maxima is a 
little higher than that obtained from taking careful account of scaling corrections arising 



from the boundary conditions, (3 c (oo) = 0.54757(63), using a Metropolis simulation 11 



The various pseudo-critical temperatures observed on differently sized lattices in a multi- 
canonical simulation fall between those observed for hot and cold starts in Metropolis 
simulations such as those performed here on the same lattices, as one might expect. 



4. Conclusions 

As Suzuki remarked many years ago [7] the anisotropic Gonihedric Ising (Fuki-Nuke) 
model may be reformulated in terms of new spin variables as a stack of decoupled 2D 
Ising models which magnetize independently. When translated back into the original 
spin variables the magnetization of each 2D Ising model is expressed in terms of a two 
spin nearest neighbour correlator (cjCTj). To take account of the independent ordering 
of planes a quasi-3-D order parameter may then be constructed by taking the modulus, 
M a b s of nearest neighbour correlators summed over each plane. 

We investigated numerically the suitability of using this and related quantities as 
an order parameter for the isotropic 3D Gonihedric Hamiltonian, H K=0 , following the 
suggestion of [6]. The transition in this case is now first order and hence no longer in 
the universality class of of the 2D Ising model, but the distinctive planar flip symmetry 
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remains. We concluded that the M^ ,z of equ. ( 19 ) might serve as an order parameter 
for the isotropic model and found that the choice of reference direction x, y, z was 
immaterial in this case. The absence of non-trivial structure in the distribution P(M2d) 
of the coplanar magnetizations suggested that such a global parameter was sufficient to 
capture the order. 

Somewhat unexpectedly, the standard magnetic susceptibility \ a l so shows the 
characteristics of an order parameter in the simulations of the isotropic Gonihedric 
model since it was measured to be one at high temperatures and dropped sharply to 
zero at the pseudo-critical point. However, \ failed to distinguish between an ensemble 
of flipped toy ground states and purely random configurations. We suggested that the 
behaviour in the full Gonihedric simulations was a dynamical consequence of freezing 
and that a multicanonical simulation might help to clarify this. 

Regarding similar order parameters in other models, the dual Gonihedric model 
may be written as strongly anisotropic Ashkin- Teller model and also possesses similar 
flip symmetries to the purely plaquette model [12]. Understanding the low temperature 
order in this dual model, whose standard magnetization and polarization are zero at all 
(3, might proceed along similar lines to those considered here. 
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